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We introduce a method for isospin restoration in the calculation of nuclear matrix elements 
(NME) for Qvjjp and decay within the framework of interacting boson model (IBM-2). With 
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I. INTRODUCTION 


The question whether neutrinos are Majorana or Dirac 
particles, and of what are their masses and phases in 
the mixing matrix remains one of the most important in 
physics today. A direct measurement of the average mass 
can be obtained from the observation of the neutrinoless 
double-/? decay (Oz^/3/3) 

—>-^±2 + 2 e ^, ( 1 ) 

two scenarios of which are shown in Fig. [T] Several ex- 



FIG. 1. Neutrinoless double-/! decay mechanism for (a) light 
neutrino exchange and (b) heavy neutrino exchange. 


Concomitant with the neutrinoless modes, there is also 
the process allowed by the standard model, 2i^/3/3, de¬ 
picted in Fig. [21 For this process, the half-life can be, to 
a good approximation, factorized in the form 



= G2v \ meC^M2i, 


( 3 ) 


The factorization here is not exact and conditions under 
which it can be done are discussed in Ref. [^. 
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FIG. 2. Double-/? decay mechanism with the emission of 2u. 
The processes depicted in Figs. [1] and |2] are of the type 
(A, Z) {A, Z -\-2) -\- 2e~ -\- anything. (4) 


periments are underway to detect this decay, and others 
are in the planning stage (for a review, see e.g. Q). The 
half-life for this decay can be written as 

[r°/ 2 ]-' = Go. iMo.l' |/(m„ , (2) 

where Go. is a phase space factor, Mq. the nuclear ma¬ 
trix element and f{mi,Uei) contains physics beyond the 
standard model through the masses rrii and mixing ma¬ 
trix elements Uei of neutrino species. 
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In recent years, interest in the processes 

{A, Z) {A, Z — 2) -\- 2e^ -h anything (5) 

has also arisen. In this case there are also the competing 
modes in which either one or two electrons are captured 
from the electron cloud (Oi^/?EC, 2i^/3EC, i?Oz/ECEC, 
2 j/ECEC). Also for these modes, the half-life can be fac¬ 
torized (either exactly or approximately) into the prod¬ 
uct of a phase space factor and a nuclear matrix element 
which then are the crucial ingredients of any double-/? 
decay calculation. 

In order to extract physics beyond the standard model, 
contained in the function / in Eq. we need an accu¬ 
rate calculation of both Go. and Mg.. These calcula¬ 
tions will serve the purpose of extracting the neutrino 
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mass (ttIj/) if OvPP is observed, and of guiding searches if 
Or^/3/3 is not observed. 

Recently we have started a systematic evaluation of 
phase space factors (PSF) and nuclear matrix elements 
(NME) for all processes of interest. The results for NME 
have beenpresented in Refs. [343 , and those for PSF in 
Refs. [3, Hill- The calculations for the nuclear matrix 
elements have been done within the framework of the 
microscopic interacting boson model (IBM-2). 

Having completed the calculations in all nuclei of in¬ 
terest, we have now readdressed them with the purpose 
of providing as accurate as possible results. As shown in 
Table XV of Ref. [3|, the Fermi matrix elements 
for 2i//3/3 decay in IBM-2 do not vanish in cases where 
protons and neutrons occupy the same major shell. Sim¬ 
ilarly, the Fermi matrix elements for Oi//3/3 decay 

are large when protons and neutrons are in the same ma¬ 
jor shell, as we can see from Table VII of Ref. where 

the quantity xf = {gv/ 9 aY/Mq^ is reported. 
The same problem with isospin was present in the quasi¬ 
particle random phase approximation (QRPA) both for 
QRPA-Tii [3| and QRPA-Jy m and it was cured re¬ 
cently Ea by changing the values of the renormalization 
constant which is adjusted in such a way as to make 

vanish. In this paper, we report on a method simi¬ 
lar in spirit, but different in practice from QRPA, and use 
it to impose the condition = 0 in IBM-2. A conse¬ 

quence of the implementation of this method is that the 
matrix elements Mp'^'^ are reduced and comparable now 
to those obtained in the Interacting Shell Model (ISM). 


i 

n,n' 

( 6 ) 

where, for s = 0, = 1, and for s = 1, = a. In 

particular, the Fermi transition operator for 2i//3/3 decay, 
obtained from Eq. ([6]) by setting A = 0, si = S 2 = 0, and 
V{r) = 1, is 

yp'"'’ (2^/3/3) (7) 

n,n' 

and for Oz^/3/3 decay obtained from Eq. ([6]) by setting 
A = 0, Si = S 2 = 0, and V{r) = H{r), is 

yp'"^ = ^ E T-+id(w) (Ojz/3^) (8) 

n,n' 

where H (r) is given in Appendix A of Ref. Q- The oper¬ 
ator o, when summed over all particles, cannot change 
isospin and its matrix elements between initial and final 
states must vanish. 

In previous IBM-2 calculations [3|, the matrix elements 
Mp’^'^ were simply discarded and the matrix elements 

Mp’''^ were kept untouched. We suggest here that a bet¬ 
ter approximation is that of modifying the mapped op¬ 
erator by imposing the condition that = 0. This 

condition can be simply implemented in our calculation 
W replacing the radial integrals of Appendix A of Ref. 
la, by 


II. FORMALISM 


The role of isospin in the Interacting Boson Model was 
extensively investigated in the 1980’s and 1990’s [T2l - [T3 | 
and is summarized in the paper "Isospin and F-spin in 
the Interacting Boson Model" by J. P Elliott [ 13 . As dis¬ 
cussed in [l3, IBM-2 wave functions have good isospin 
in heavy nuclei with a neutron excess. The problem 
arises only in light nuclei where protons and neutrons 
occupy the same orbits {sd and pf shells). For these nu¬ 
clei one needs to introduce an isospin invariant form of 
IBM, called IBM-3 (T3|. IBM-2 can still be used in light 
nuclei if the parameters of the interaction are obtained 
by projection from those of the isospin invariant IBM-3. 
(For the nuclei discussed in this paper, only ^®Ca and 
^®Ti are in a situation in which IBM-3 or a projected 
form should be used.) As a result, isospin does not pose 
a problem for IBM-2 wave functions. 

The problem with isospin in IBM-2 arises in the map¬ 
ping of the fermion operator for Oz//3/3 and 2z^/3/3 decay. 
The expression for the fermionic transition operator of 
type Fermi (F), Gamow-Teller (GT), and tensor (T) is 


2i//3/? : 

^ki ,0^k2,0^k,0^X,0^ji ,j'^ ^j2 J 2 ^n2,n2 ^l2 ,^2 


0z//3^ : 

,0^/c2 ,0^A:,0^A,0^_7 i j'J ^j2 ,^2 ^n2 ^^ 2 X 2 

X i?^°’°’°^(ni,Zi,n2,Z2,'rei,/i,n2,Z2) 


( 10 ) 


where 


Rff'^\ni,li,n2, /2, n'l, Zj, 712, ^ 2 ) 


1 


/o Pp{p + A) 


p^dp 


sin(pri) 2 , 


X / Rnih {'ri)Rn'^v^ {ri) ' '2 "' rfdri 

Jo Pi 

X [ Rn-,iXr2)Rn'^V^{r2f-^^^^^^^rldr2. 

Jo P 2 

This prescription will guarantee that the F matrix el- 
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ements vanish for 2i//3/3, and will reduce the F matrix 
elements for Qv(3j3 by subtraction of , which is the 

monopole term in the expansion of the matrix element 
into multipoles. It is similar to the prescription used in 
Ref. (see Fig. 4 of Ref. M)- 

Although the method described in this section is not 
an isospin restoration of the IBM-2 wave functions which 
have already good isospin but rather a restoration of the 
isospin properties of the mapping of the transition oper¬ 
ator, we shall, nonetheless, for simplicity refer to it in the 
following sections as "simply isospin restoration". 

III. RESULTS FOR OuPP 

From here on, the calculation of the matrix elements 
in the interacting boson model proceeds in the same way 
as in Refs. We do not repeat the formulas but only 

report the results of the calculation. In the calculation 
one needs to specify the parametrization of short-range 
correlations (SRC). In earlier calculations Q the Miller- 
Spencer parametrization was used. It has now become 
clear that Argonne and CD-Bonn parametrizations are 
more appropriate. Here we use throughout the Argonne 
parametrization of the correlation function 

/(r) = l-ce“’''(I-&r2), (12) 

that is a = 1.59fm“^, b = I.45fm“^, and c = 0.92. We 
write 

Mo. = 

+ Mp’'\ 

with the ratio gy / Qa explicitly displayed in front of 



A. OvPP decay with light neutrino exchange 

In Table HI we show the results of our calculation of 
the nuclear matrix elements to the ground state, 0)^, and 
the first excited state, Oj, broken down into GT, F, and 
T contributions and their sum according to Eq. (HU. 
The parameters of the IBM-2 Hamiltonian used in this 
calculation are those given in Table XXHI of Ref. @ 
(with the exception of ^®'^Gd whose parameters are given 
in Ref. [Hi). 

When cornpared with the matrix elements without the 
restoration [ 3 , we see a considerable reduction of the 
F matrix elements to values comparable to those of the 
shell model. The overall reduction in is ~ 15%. 

Our results are compared with QRPA-Tii with isospin 
restoration (Argonne SRC) [IH and ISM (UCOM SRC) 
[H | in Table ITTl and Fig.[3l 



FIG. 3. (Color online) IBM-2 isospin restored results for 
0i^~/3~ decay compared with QRPA-Tii flll| and the ISM 

The reduction in the Fermi matrix elements 
brought in by the isospin restoration is shown in Ta¬ 
ble Hill where the quantity xf = {gv/ 9 aYMqP 
is shown for the old and new calculation and compared 
with QRPA without (old )and with (new) isospin restora¬ 
tion, and with ISM. Our isospin restored Fermi matrix 
elements are comparable to those of the ISM, but a fac¬ 
tor of 2-3 smaller than the isospin restored QRPA-Tii 
results. This may be due to the fact that both in IBM-2 
and ISM the model space is rather restricted, while in 
QRPA several major shells are included. 


B. OuPP decay with heavy neutrino exchange 

These matrix elements can be simply calculated by 
replacing the potential v{p) = 2Tr~^[p{p + A)\~^ in 
jj(fci,fc 2 ,A) potential Vh{p) = 27r“^(memp)“^. Ta¬ 

ble HY] gives the corresponding matrix elements. The in¬ 
dex "h" distinguishes these matrix elements from those 
with light neutrino exchange. Our results are compared 
with results of QRPA-Tii [2l| and ISM [H] in Table Ivl 


1. Sensitivity to parameter ehanges, model assumptions 
and operator assumptions 

The sensitivity of IBM-2 to parameter changes, model 
assumptions, and o^rator assumptions was discussed in 
great detail in Ref. 0. We do not repeat this discussion, 
but only note that because of isospin restoration the sen¬ 
sitivity of F matrix elements to isospin purity decreases 
from 10 % to 1 %, including the case of ^®Ca decay, which 
previously was treated separately from the others. Our 
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TABLE 1. IBM-2 nuclear matrix elements (dimensionless) for QvP /3 decay with Argonne SRC, Qv/dA ~ 1/1.269, and 

isospin restoration. 


A 



0 ^ 




0 ^ 










48Ca 

1.73 

-0.30 

-0.17 

1.75 

3.78 

-0.27 

-0.12 

3.82 

76Ge 

4.49 

-0.68 

-0.23 

4.68 

1.95 

-0.27 

-0.09 

2.02 

82Se 

3.59 

-0.60 

-0.23 

3.73 

0.92 

-0.13 

-0.05 

0.95 

96 Zr 

2.51 

-0.33 

0.11 

2.83 

0.04 

-0.01 

0.00 

0.05 

looMo 

3.73 

-0.48 

0.19 

4.22 

0.99 

-0.13 

0.05 

1.12 

llOpd 

3.59 

-0.40 

0.21 

4.05 

0.46 

-0.05 

0.03 

0.52 

ii^Cd 

2.76 

-0.33 

0.14 

3.10 

0.84 

-0.09 

0.03 

0.93 

i24Sn 

2.96 

-0.57 

-0.12 

3.19 

2.21 

-0.41 

-0.09 

2.38 

128Xg 

3.80 

-0.72 

-0.15 

4.10 

2.65 

-0.47 

-0.09 

2.85 

laopg 

3.43 

-0.65 

-0.13 

3.70 

2.52 

-0.45 

-0.08 

2.71 

i34Xe 

3.77 

-0.68 

-0.15 

4.05 

2.19 

-0.36 

-0.06 

2.35 

i36Xe 

2.83 

-0.52 

-0.10 

3.05 

1.49 

-0.24 

-0.03 

1.60 

148Nd 

2.00 

-0.38 

0.07 

2.31 

0.25 

-0.05 

0.01 

0.29 

150Nd 

2.33 

-0.39 

0.10 

2.67 

0.40 

-0.06 

0.02 

0.45 

i54Sm 

2.49 

-0.36 

0.11 

2.82 

0.37 

-0.04 

0.01 

0.41 

160 Gd 

3.64 

-0.45 

0.17 

4.08 

0.76 

-0.11 

0.04 

0.87 

198pt 

1.90 

-0.33 

0.09 

2.19 

0.08 

-0.02 

0.01 

0.10 

232Th 

3.58 

-0.44 

0.18 

4.04 

0.12 

-0.02 

0.01 

0.15 

238^ 

4.27 

-0.53 

0.21 

4.81 

0.34 

-0.05 

0.02 

0.40 


error estimate for Qvjpp is now 16% for all nuclei. In the 
case of OvhPP we also estimate a reduced sensitivity of F 
matrix elements from 10% to 1%, and a reduced sensitiv¬ 
ity of the matrix elements F + GT to SRC from 50% to 
25%. This sensitivity is estimated by comparing matrix 
elements with Argonne-CD-Bonn and UCOM SRC. Our 
total error estimate for OvhPfi is now 28% for all nuclei. 
Our final matrix elements, with error estimate are given 
in Table IvH 

In addition to IBM-2, QRPA, and ISM, other cal¬ 
culations have been recently done. In Fig. 2] we com¬ 
pare our results with all available calculations including 
density functional theory (DFT) [2^ and Hartree-Fock- 
Bogoliubov (HFB) theory (^ . 


IV. RESULTS FOR QvP+P+ 

Matrix elements for double positron (/3+/3+) decay and 
the related processes {EC(3'^) and {ECEC) can be cal¬ 
culated in a similar way. 


A. QvP ^and related processes with light 
neutrino exchange 

In Table IVIII we show the results of our calculation of 
the matrix elements to the ground state, 0)^, and first 
excited, 0^, state broken down into GT, F, and T con¬ 
tributions and their sum according to Eq. (USD. 



FIG. 4. (Color online) IBM-2 (Argonne) results for QvP P 
nuclear matrix elements compared with QRPA-Tii (Argonne) 
[H, ISM (UCOM) [l^QRPA-Jy (UCOMlil,Ugl, QRPA- 
deformed (CD-Bonn)l23|, DFT (UCOM) gjand HFB (M- 

s) m. 


The parameters of the IBM-2 Hamiltonian used in this 
calculation are those in Tables II and VI of Ref. in, 
respectively. Also here, as in the previous Sect. Ill, we 
see that the F matrix elements are considerably reduced 
by isospin restoration in comparison with those without 
restoration given in Table VIII of Ref. This is also 
seen in Table IVIIII where the quantity xp is shown. 






















5 


TABLE 11. Comparison among nuclear matrix elements for 
OuP~P~ decay to ground state, 0]*', in IBM-2 with Argonne 
SRC, QA = 1.269, and isospin restoration, QRPA-Tii with Ar¬ 
gonne SRC, qa = 1.27, and isospin restoration [H, and ISM 
with UCOM SRC and qa = 1-25 [2^. All matrix elements 
are in dimensionless units. 


Decay 




IBM-2 

QRPA-Tii 

ISM 

48Ca^48Tj 

1.75 

0.54 

0.85 

76Ge^76gg 

4.68 

5.16 

2.81 

82Se^®2Kr 

3.73 

4.64 

2.64 

96Zr^-9®Mo 

2.83 

2.72 


io°Mo^i°ORu 

4.22 

5.40 


iiopd^iiocd 

4.05 

5.76 



3.10 

4.04 


124Sn^l24Te 

3.19 

2.56 

2.62 

i28Te^.i28xe 

4.10 

4.56 

2.88 

i30Te^.i30xe 

3.70 

3.89 

2.65 

i34Xe^-i34Ba 

4.05 



i36Xe^-i36Ba 

3.05 

2.18 

2.19 

148Nd^l48Sm 

2.31 



isoNd^isogm 

2.67 



i54Sm^-i54Gd 

2.82 



i60Gd^-i60Dy 

4.08 



198pt^l98jjg 

2.19 



232Th_ 

4.04 



238u^238p^ 

4.81 




TABLE III. Comparison between Fermi matrix eleme nts, yf , 
for 0uP~P~ decay in IBM-2, QRPA-Tii [Till and ISM [ig. l^ . 


Xf 


Decay 

IBM-2 

QRPA-Tii 

ISM 


old 

new 

old 

new 


48Ca 

-0.39 

-0.11 

-0.93 

-0.32 

-0.14 

™Ge 

-0.37 

-0.09 

-0.34 

-0.21 

-0.10 

82 Se 

-0.40 

-0.10 

-0.35 

-0.23 

-0.10 

96zr 

-0.08 

-0.08 

-0.38 

-0.23 


looMo 

-0.08 

-0.08 

-0.30 

-0.30 


llOpd 

-0.07 

-0.07 

-0.33 

-0.27 

-0.16 

ii^Gd 

-0.07 

-0.07 

-0.30 

-0.30 

-0.19 

i24Sn 

-0.34 

-0.12 

-0.40 

-0.27 

-0.13 

i28Te 

-0.33 

-0.12 

-0.38 

-0.27 

-0.13 

i30Te 

-0.33 

-0.12 

-0.39 

-0.27 

-0.13 

i34Xe 


-0.11 




i36Xe 

-0.11 

-0.11 

-0.38 

-0.25 

-0.13 

148Nd 

-0.12 

-0.12 




150Nd 

-0.10 

-0.10 




i54Sm 

-0.09 

-0.09 




i™Gd 

-0.08 

-0.08 




198pt 

-0.11 

-0.11 




232Th 

-0.08 

-0.08 




238^ 

-0.08 

-0.08 





V. RESULTS FOR 2vPP 


Our results are compared with other available calcula¬ 
tions in Table IdO For /3+/3+, ECP^, and EC EC decay 
there are no QRPA calculations with isospin restoration 
and thus the comparison is only meant to show the re¬ 
duction in the F matrix element in IBM-2 brought in by 
isospin restoration. 


B. OvP^p^ and related processes with heavy 
neutrino exchange 

These matrix elements are obtained in the same way 
as in Sect. III.2 and are given in Table DO 


1. Sensitivity to parameter ehanges, model assumptions 
and operator assumptions 


Isospin restoration has a major consequence on matrix 
elements for 2vPP decay, since F matrix elements van¬ 
ish when isospin restoration is imposed. 2z//3/3 matrix 
elements can be easily calculated in IBM-2 using closure 
approximation (CA). In this approximation the matrix 
elements M 2 v^ which appear in the half-life Eq. m can 
be written as 


where 




M^gt = ( 0^ 


\m^gt 

( 

1- 

cs 

Agt 

VSA/ 

Ap 


■ Or,' 




= I 0^ 


E 


t t 

T T / 
'n' n' 


0} 


(14) 


(15) 


The closure energies Agt and Ap are defined by 


The sensitivity here is identical to that described in 
Sect. Ill for 0i'P~P~. Our final matrix elements with 
error estimate are given in Table IXII 


Agt = 2 {QfiP + 2TOeC^) + (Ei+^n) — Ep 
Ap = — {Qpp + 2meC^) + (E-o+.Af) ~ Ep 


(16) 
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TABLE IV. Nuclear matrix elements for the heavy neutrino exchange mode of the neutrinoless double-/?” decay to the ground 
state (columns 2,3,4) and 5) and to the first excited state (columns 6,7,8, and 9) using the microscopic interacting boson model 
(IBM-2) with isospin restoration and Argonne SRC and gvjgA = 1/1.269. 


A 


0^ 



0^ 










48Ca 

53.5 

-23.2 

-21.3 

46.6 

44.8 

-8.8 

- 6.5 

43.7 

76Ge 

104 

-42.8 

-26.9 

104 

38.6 

-14.9 

-9.8 

38.1 

82Se 

87.2 

-37.1 

-27.3 

82.9 

16.8 

-6.5 

-4.6 

16.2 

96 Zr 

67.9 

-29.2 

12.7 

98.7 

1.4 

-0.6 

0.3 

2.1 

looMo 

111 

-46.8 

24.2 

164 

29.3 

-12.4 

6.4 

43.3 

llOpd 

100 

-41.4 

27.7 

154 

13.5 

-5.6 

3.8 

20.9 

iiecd 

73.9 

-31.2 

16.9 

no 

18.0 

-7.5 

3.5 

26.1 

i24Sn 

73.7 

-33.1 

-14.9 

79.3 

50.1 

-22.2 

-9.9 

54.0 

i28xe 

93.4 

-41.7 

-18.3 

101 

55.7 

-24.4 

-10.3 

60.5 

i30xe 

84.8 

-37.9 

-16.6 

91.8 

51.5 

-22.6 

-9.3 

56.2 

i34Xe 

86.6 

-39.3 

-19.8 

91.2 

38.7 

-17.3 

-7.9 

41.5 

i36Xe 

66.8 

-29.7 

-12.7 

72.6 

25.6 

-11.0 

-4.1 

28.3 

148 Nd 

72.8 

-32.7 

9.6 

103 

8.1 

-3.7 

1.0 

11.4 

isoNd 

81.1 

-35.6 

13.2 

116 

12.2 

-5.3 

1.8 

17.3 

i54Sm 

78.1 

-33.7 

13.8 

113 

8.9 

-3.8 

1.2 

12.4 

160 Gd 

106 

-44.6 

21.5 

155 

26.7 

-11.4 

6.2 

40.0 

198pt 

71.4 

-31.9 

12.8 

104 

4.0 

-1.8 

0.9 

6.1 

232Th 

107 

-44.0 

24.4 

159 

6.2 

-2.7 

1.9 

9.9 

238^ 

127 

-52.5 

28.7 

189 

12.7 

-5.4 

3.4 

19.4 


where is a suitably chosen excitation energy in 

the intermediate odd-odd nucleus. The matrix elements 
can be simply calculated by replacing the neutrino 
potential v{p) 


V2u{p) 


Hp) 

n2 


(17) 


which is the Fourier-Bessel transform of the configuration 
space potential V{r) = 1. 

In order to confirm that in isospin-restored IBM-2 cal¬ 
culation the Fermi matrix elements for 2vj3j3 decay van¬ 
ish we have calculated and The results are 

given in Table Ixlll We can see from this table that 
indeed vanish. The small values ^ 0.01 are an indication 
of our numerical accuracy in calculating the radial inte¬ 
grals 

Using the results in Table IXlIl one can redo the analysis 
of Ref. 0 and extract the values of the effective gA,eff 
from 


^e// ^ / 9A^\ 

\ 9A J 

with extracted from experiment (Slj as compiled 

in Ref. [ 3 . The corresponding results are shown in Fig. [S] 
Isospin restoration has no effect on the extracted values of 
gA,eff, since in the previous analysis @ the Fermi matrix 
elements were simply discarded. The difference 


between Fig. 5 of this paper and Fig. 13 of is only 
due to the fact that we have used Argonne SRC instead 
of Miller-Spencer. 



Mass number 


FIG. 5. (Color online) Value of gA,eff extracted from exper¬ 
iment for IBM-2 and ISM. 


We also note that, very recently, gA,eff values have 
also been extracted in QRPA-Tii and QRPA-Jy 
with results similar to those in Fig. O 
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TABLE V. Neutrinoless double-/?” decay nuclear matrix ele¬ 
ments to ground state, 0/*", in IBM-2 with isospin restoration, 
Argonne SRC and gv j9 a = 1/1.269, QRPA-Tii with isospin 
restoration and Argonne SRC and ISM with UCOM SRC 
for the heavy neutrino exchange mode. All matrix ele¬ 
ments are in dimensionless units. 


Decay 




IBM-2 

QRPA-Tii 

ISM 

48Ca^48Tj 

46.6 

40 

47.5 

76Ge^76gg 

104 

287 

138 


82.9 

262 

127 

96Zr^.96Mo 

98.7 

184 


iooMo^-1°°Ru 

164 

342 


iiopd^iiocd 

154 

333 



no 

209 


124gj^^l24pg 

79.3 

184 


i28Te^.i28xe 

101 

302 


i30Te^.i30xe 

91.8 

264 


i34Xe^-i34Ba 

91.2 




72.6 

152 


148Nd^l48gj^ 

103 



i50Nd^i59Sm 

116 




113 



leoQd^isoDy 

155 



198pt^l98jjg 

104 



232Th_ 

159 



238u^238p^ 

189 




VI. RESULTS FOR 2iy/3+/3+ AND COMPETING 
MODES 


Matrix elements for 2i//3+/3+ and related processes can 
be obtained in the same way as in Sect. V. The results 
are given in Table IXIIII 


VII. EXPECTED HALF-LIVES AND LIMITS ON 
NEUTRINO MASS 


A. Light neutrino exchange 


The calculation of nuclear matrix elements in IBM-2 
with isospin restoration can now be combined with phase- 
space factors sai to produce our final results for half- 
lives for light neutrino exchange in Table IXTVl and Fig.|6l 

For light neutrino exchange 


\f{mi,Uei)f 


{rriu) 

me 


(19) 


The average light neutrino mass is constrained by atmo¬ 
spheric, solar, reactor and accelerator neutrino oscillation 


TABLE VI. Final double-/? decay IBM-2 matrix elements 
with isospin restoration, Argonne SRC and error estimate. 

Decay Light neutrino exchange Heavy neutrino exchange 


48Ca 

1.75(28) 

47(13) 

76Ge 

4.68(75) 

104(29) 

82 Se 

3.73(60) 

83(23) 

96zr 

2.83(45) 

99(28) 

looMo 

4.22(68) 

164(46) 

llOpd 

4.05(65) 

154(43) 

“6cd 

3.10(50) 

110(31) 

i24Sn 

3.19(51) 

79(22) 

i28Te 

4.10(66) 

101(28) 

i30Te 

3.70(59) 

92(26) 

i34Xe 

4.05(65) 

91(26) 

i36Xe 

3.05(59) 

73(20) 

148Nd 

2.31(37) 

103(29) 

150Nd 

2.67(43) 

116(32) 

i84Sm 

2.82(45) 

113(32) 

leoQd 

4.08(65) 

155(43) 

198pt 

2.19(35) 

104(29) 

232Th 

4.04(65) 

159(45) 

238^ 

4.81(77) 

189(53) 



FIG. 6. (Color online) Expected half-lives for {rtiv) = 1 eV, 
QA = 1.269 and IBM-2 isospin restored nuclear matrix ele¬ 
ments. The points for ^^®Te, ^®^Xe, and ^'^^Nd decays are not 
included in this figure. The figure is in semilogarithmic scale. 


experiments to be 

(m^) = |ci3Ci2mi + Si3m3e"‘^^| , 


Cij cos '0 ij , 
( 2 2 2 \ 


Sij 

m? + mo / 


<P2,3 = [0,27r], 

5m? 5m? 

— 


±Am' 


( 20 ) 


























TABLE VII. Nuclear matrix elements (dimensionless) for neutrinoless ECP'^ , and EC EC decays with Argonne 

SRC and gv/Q a = 1/1.269, in IBM-2 with isospin restoration. 


Nucleus 



0 ^ 




0 ^ 


M^ot 








58Ni 

2.33 

-0.23 

0.15 

2.61 

2.21 

- 0.20 

0.10 

2.44 

®4Zn 

5.22 

-0.61 

-0.16 

5.44 

0.68 

-0.06 

- 0.02 

0.70 

78Kr 

3.79 

-0.61 

-0.24 

3.92 

0.87 

-0.14 

-0.06 

0.90 

96Ru 

2.51 

-0.37 

0.11 

2.85 

0.03 

- 0.01 

0.00 

0.04 

loecd 

3.16 

-0.38 

0.19 

3.59 

1.55 

-0.16 

0.08 

1.72 

i24Xe 

4.42 

-0.82 

-0.19 

4.74 

0.74 

-0.14 

-0.03 

0.80 

i30Ba 

4.36 

-0.80 

-0.18 

4.67 

0.32 

-0.06 

- 0.01 

0.34 

i36Ce 

4.23 

-0.76 

-0.16 

4.54 

0.35 

-0.06 

- 0.01 

0.38 

156Dy 

2.80 

-0.40 

0.13 

3.17 

1.53 

-0.23 

0.08 

1.75 

i64Er 

3.46 

-0.44 

0.22 

3.95 

1.02 

- 0.10 

0.05 

1.13 


4.12 

-0.57 

0.20 

4.67 

0.26 

-0.05 

0.02 

0.31 


TABLE VIII. Ratio Fermi to Gamow-Teller matrix elements, 
Xf, for neutrinoless , ECP^ , and ECEC in the IBM-2 
with isospin restoration compared with available QRPA re¬ 
sults. 




XF 


Decay 


IBM-2 

QRPA^ 


old 

new 


58Ni 

-0.06 

-0.06 

-0.14 

64 Zn 

-0.31 

-0.07 


78Kr 

-0.38 

- 0.10 

-0.27 

66 Ru 

-0.09 

-0.09 

-0.23 

loscd 

-0.07 

-0.07 

-0.23 

i24Xe 

-0.34 

- 0.12 

-0.23 

i30Ba 

-0.32 

- 0.11 

-0.23 

436Ce 

-0.32 

- 0.11 

-0.26 

156Dy 


-0.09 


i64Er 


-0.08 


180-W 


-0.09 



Reference l2Sll . No isospin restoration. 


TABLE IX. IBM-2 matrix elements with Argonne SRC and 
isospin restoration for neutrinoless P'^P'^, ECP'^ , and ECEC 
compared with available QRPA calculations. 


Decay 


0 ^ 

oj 

IBM-2 

QRPA^ 

IBM-2 QRPA 

58Ni 

2.61 

1.55 

2.44 

64 Zn 

5.44 


0.70 

78Kr 

3.92 

4.16 

0.90 

66Ru 

2.85 

3.23 4.29^ 0.04 2.31^ 

106Cd 

3.59 

4.10 7.54^= 1.72 0.6U 

I24xe 

4.74 

4.76 

0.80 

i30Ba 

4.67 

4.95 

0.34 

i36Ce 

4.54 

3.7 

0.38 

156Dy 

3.17 


1.75 

i64Er 

3.95 


1.13 

iso^Y 

4.67 


0.31 


Reference 
^ Reference 
Reference 


g . No isospin restoration. 

(UCOM SRC). No isospin restoration. 
(UCOM SRC). No isospin restoration. 


Using the best fit values [s^l 

sin^ Pi2 = 0.213, sin^ = 0.016, 
sin^ P 23 = 0.466, Sm^ = 7.67 x 10"® eV^ (21) 
Am^ = 2.39 X 10"^ eV^ 

we obtain the values given in Fig. [T] In this figure we 
have added the current limits, for = 1.269, com¬ 
ing from CUORICINO |3^ I GEX 0, NEMO-3 [13, 
KamLAND-Zen [H, EXO [Hi and GERDA 0 exper¬ 
iments. Also, henceforth we use c=l to conform with 
standard notation. 


B. Heavy neutrino exchange 

The half-lives for this case are calculated using the for¬ 
mula 

— / -i\ rrr (22) 

V = mp{m,^)= (C/efeJ —. 

k— heavy ^ 

The expected half-lives for \r]\ = 10“^, and using the 
IBM-2 matrix elements of Table are shown in Ta¬ 
ble [XV] For other values of 77 they scale as \ri\‘^. There 
are no direct experimental bounds on 77 . Recently, Tello 
et al. [13 have argued that from lepton flavor violating 
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TABLE X. Nuclear matrix elements (dimensionless) for heavy neutrino exchange for neutrinoless 13'^ / EC / EC EC decay 

in IBM-2 with isospin restoration, Argonne SRC, and gv/Q a = 1/1.269. 


Nucleus 



0 ^ 




0 ^ 










58Ni 

55.1 

-23.1 

18.6 

88.0 

36.3 

-15.8 

8.33 

54.5 

64 Zn 

103 

-38.9 

-18.5 

109 

lO.l 

-3.20 

- 2.00 

10.1 

78Kr 

89.8 

-38.5 

-30.6 

83.1 

21.1 

-9.12 

-7.22 

19.5 

96Ru 

67.5 

-30.6 

12.5 

99.0 

0.32 

-0.08 

0.32 

0.59 

106Cd 

87.8 

-38.1 

26.5 

138 

34.0 

-14.7 

8.75 

51.9 

I24xe 

105 

-47.9 

-25.0 

no 

18.1 

-8.24 

-4.31 

18.9 

i30Ba 

103 

-46.4 

-23.7 

108 

8.07 

-3.68 

-1.90 

8.45 

436Ce 

95.8 

-43.2 

- 21.8 

101 

8.24 

-3.73. 

-1.89 

8.66 

156Dy 

82.6 

-37.0 

17.5 

123 

47.6 

-21.4 

10.4 

71.3 

i64Er 

108 

-46.8 

32.9 

170 

23.6 

-9.95 

5.96 

35.8 

180w 

119 

-53.3 

28.1 

180 

10.7 

-4.85 

2.91 

16.6 


TABLE XL Final ECj3+, and EC EC IBM-2 matrix 

elements with isospin restoration, Argonne SRC, and their 
error estimate. 


Decay Light neutrino exchange Heavy neutrino exchange 


58Ni 

2.61(42) 

88(25) 

64 Zn 

5.44(87) 

109(31) 

78Kr 

3.92(63) 

83(23) 

66 Ru 

2.85(46) 

99(28) 

106Cd 

3.59(57) 

138(39) 

I24xe 

4.74(76) 

110(31) 

i30Ba 

4.67(75) 

108(30) 

436Ce 

4.54(73) 

101(28) 

156 J3y 

3.17(51) 

123(34) 

i64Er 

3.95(63) 

170(48) 

isow 

4.67(75) 

180(50) 


processes and from large hadron collider (LHC) exper¬ 
iments one can put some bounds on the right-handed 
leptonic mixing matrix Ugk,heavy and thus on rj. In the 
model of Ref. |43| . when converted to our notation, r] can 
be written as 


V = 


^WR 


k—heavy 


rup 
mkh ’ 


(23) 


where Mw is the mass of the LL-boson, Mw = (80.41 ± 
0.10) GeV [3, Mwr is the mass of IRi?-boson, and 
V = {M\yr/M\y)‘^ U. By comparing the calculated half- 
lives with their current experimental limits, we can set 
limits on the lepton nonconserving parameter I 77 I, shown 
in Table IXVl 

If we write 


V = 


^WR 


TABLE XII. 2vf3~P~ matrix elements (dimensionless) to the 
ground state (columns 2 and 3) and to the first excited state 
(columns 4 and 5) using the microscopic interacting boson 
model (IBM-2) with isospin restoration and and Argonne SRC 
in the closure approximation. 


Nucleus 


0 ^- 


0 ^ 

M^ar^ 


M^^'> 


48Ca 

1.64 

- 0.01 

5.07 

- 0.01 

™Ge 

4.44 

- 0.01 

2.02 

- 0.00 

82Se 

3.59 

- 0.01 

1.05 

- 0.00 

96zr 

2.28 

- 0.00 

0.04 

- 0.00 

looMo 

3.05 

- 0.00 

0.81 

- 0.00 

llOpd 

3.08 

- 0.00 

0.38 

- 0.00 

“ 6 Cd 

2.38 

- 0.00 

0.83 

- 0.00 

i24Sn 

2.86 

- 0.01 

2.19 

- 0.00 

i28Te 

3.71 

- 0.01 

2.70 

- 0.00 

i30Te 

3.39 

- 0.01 

2.64 

- 0.00 

434Xe 

3.69 

- 0.01 

2.34 

- 0.00 

i36Xe 

2.82 

- 0.01 

1.65 

- 0.00 

148Nd 

1.31 

- 0.00 

0.18 

- 0.00 

150Nd 

1.61 

- 0.00 

0.31 

- 0.00 

464Sm 

1.95 

- 0.00 

0.35 

- 0.00 

160Gd 

3.08 

- 0.00 

0.53 

- 0.00 

198pt 

1.06 

- 0.00 

0.03 

- 0.00 

232Th 

2.75 

- 0.00 

0.08 

- 0.00 

23S]j 

3.35 

- 0.00 

0.24 

- 0.00 


we can also set limits on the average heavy neutrino mass, 
). This limit is model dependent. In Ref. [d^ a value 
of Mwr = 3.5 TeV was used. We use here instead a lower 
value Mwr = 1.75 TeV, obtaining the limits on (rrii,^) 
shown in the last column of Table IXVl For other values 
of Mwr it scales as 


(24) 
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TABLE XIII. 2ul3+l3+, 2uECI3+, and 2vECEC nuclear ma¬ 
trix elements (dimensionless) to the ground state (columns 2 
and 3) and to the first excited state (columns 4 and 5) using 
IBM-2 with isospin restoration and and Argonne SRC in the 
closure approximation. 


Nucleus 


0^ 


0^ 



M'ff 


58Ni 

2.11 

-0.00 

2.34 

-0.00 

64 Zn 

5.20 

-0.01 

0.71 

-0.00 

78Kr 

3.67 

-0.01 

0.83 

-0.00 

96Ru 

2.17 

-0.00 

0.05 

-0.00 

106Cd 

2.57 

-0.00 

1.47 

-0.00 

424Xe 

4.24 

-0.01 

0.71 

-0.00 

430Ba 

4.22 

-0.01 

0.30 

-0.00 

i36Ce 

4.17 

-0.01 

0.34 

-0.00 

156Dy 

2.20 

-0.00 

1.15 

-0.00 

i64Er 

2.58 

-0.00 

0.90 

-0.00 

180 w 

3.09 

-0.01 

0.12 

-0.00 



FIG. 7. (Color online) Current limits to {rriv) from CUORI- 
CINO Jp, IGEX Hi, NEMO-3 [13, KamLAND-Zen [H, 
EXO I39II . and CERDA [iOl , and IBM-2 Argonne SRC isospin 
restored nuclear matrix elements and qa = 1.269. The value 
of Ref. ( 4 ^ is shown by X. The figure is in logarithmic scale. 


TABLE XIV. Left: Calculated half-lives in IBM-2 Argonne 
SRC for neutrinoless double-/! decay for {rriv) = 1 eV and 
qa = 1.269. Right: Upper limit on neutrino mass from 
current experimental limit from a compilation of Barabash 

a The value reported W Klapdor-Kleingrothaus et al. 

, IGEX collaboration 13^. and the recent limits from 
KamLAND-Zen fs^ . EXO Isj . and CERDA [iOl are also 
included. 


Decay 

rO;,( 1024 yr) 

^D 2 .exp(yr) 

{m^) (eV) 

48 Ca^ 48 Ti 

1.33 

> 5.8 X 10^2 

< 4.8 

™Ge^™Se 

1.95 

> 1.9 X 1026 

< 0.32 



1.2 X 1026 _^ 

0.40 



> 1.6 X 1026 ^ 

< 0.35 



> 2.1 X 1026 <= 

< 0.30 

82 Se^® 2 Er 

0.71 

> 3.6 X 1023 

< 1.4 

96 Zr^ 96 ]y[Q 

0.61 

> 9.2 X 1021 

< 8.1 

406 Mo^ 4 ™Ru 

0.36 

> 1.1 X 1024 

< 0.57 

iiopd_^iiocd 

1.27 



446 Cd-)> 446 Sn 

0.63 

> 1.7 X 1023 

< 1.9 

124 Sn^l 24 Tg 

1.09 



428 Te-)> 428 Xe 

10.19 

> 1.5 X 1024 

< 2.6 

i 30 Te_^i 30 xe 

0.52 

> 2.8 X 1024 

< 0.43 

i 34 xe^i 24 Ba 

10.23 



436 Xe^ 436 Bj^ 

0.74 

> 1.9 X 102 _^ 

< 0.20 



> 1.6 X 1026 « 

< 0.22 

148 xd-^ 148 Sni 

1.87 



150 Nd-). 150 Sjn 

0.22 

> 1.8 X 1022 

< 3.5 

464 Sm^ 464 Q^ 

4.19 



460 Gd^ 460 j 3 y 

0.63 



198 pt_^ 198 Hg 

2.77 



232 Th_^ 232 u 

0.44 



238 u^ 238 p^ 

0.13 




=>• Ref. 
Ref. 
Ref. 
Ref. 
® Ref. 




C. Expected half-lives for double positron decay 

Although no limits are available here, we include for 
completeness in Fig. [8] our expected half-lives for double 
positron decay and positron emitting electron capture. 

The matrix elements reported in Ref. [3| for 
ROvECEC were already obtained with isospin restora¬ 
tion and therefore expected half-lives for this process are 
not reported here. 


If both light and heavy neutrino exchange contribute, 
the half-lives are given by 



— C' 


(0) 

Ou 



{m„) 

me 


2 




(25) 


D. Effect of renormalization of qa on expected 
half-lives 

Half-lives of decay depend on qa as g\. In Fig. 

5 we have shown the values of gA,eff both for IBM-2 
with isospin restoration and for ISM as extracted from 
2v(3l3 decay. There is much discussion at the present time 
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TABLE XV. Left: Calculated half-lives for neutrinoless double-/! decay with exchange of heavy neutrinos for rj = 2.75 x 10“^ and 
QA = 1.269. Right: Upper limits of I 77 I and lower limits of heavy neutrino mass (see text for details) from current experimental 
limit from a compilation of Barabrish 0- The value reported by Klapdor-Kleingrothaus et al. [421 . IGEX collaboration (3^ . 
and the recent limit from KamLAND-Zen [s^l. EXO (3^. and GERDA (iOl are also included. 


Decay 



|r/|(10-3) 

(m.J(GeV) 

4®Ca^-'‘®Ti 

0.72 

> 5.8 X 10^2 

<0.36 

>11.9 

76Ge^76se 

1.51 

> 1.9 X lO^s 

<0.028 

>148 



1.2 X 1025_^ 

0.035 

118 



> 1.6 X 102’(^ 

<0.031 

>136 



> 2.1 X 1028<= 

< 0.027 

156 

82Se^82Kr 

0.55 

> 3.6 X 1023 

<0.12 

>34 

96Zr^96Mo 

0.19 

> 9.2 X 1021 

<0.46 

>9.15 

iooMo->i™Ru 

0.09 

> 1.1 X 1024 

<0.028 

>146 

iiopd-^iiocd 

0.33 





0.19 

> 1.7 X 1023 

<0.11 

>39.5 

124Sn^l24Tg 

0.67 




i28Te->i28xe 

6.43 

> 1.5 X 1024 

<0.21 

>20.2 

i30Te_^i30xe 

0.32 

> 2.8 X 1024 

<0.034 

>123 


8.57 





0.50 

> 1.9 X 1023_f_ 

<0.016 

>257 



> 1.6 X 1023® 

<0.018 

>236 

148Nd^l48Sni 

0.36 




isoNd^isogjj^ 

0.05 

> 1.8 X 1022 

<0.16 

>26.3 

454Sm^-454Gd 

1.00 




460Gd^460j3y 

0.17 




198pt^l98Hg 

0.48 




232Th->.232u 

0.11 




238U^238p^ 

0.03 






on whether or not 0i//?/3 is equally quenched as 2i//3/3 
(and single /3-decay). In order to investigate the possi¬ 
ble impact of quenching of gA we present in Table IXVIl 
the predicted half-lives under the assumption of maximal 
quenching 


= 1.269A-°18, 

= 1.269A-°'^2 


(26) 


erator in the calculation of the Fermi matrix elements 
within of the framework of IBM-2, and done a consis¬ 
tent calculation of Oi^/3/3, OvhPP, and 2z//3/3 NME in the 
closure approximation. With this method, the F matrix 
elements for 2i//3/3 decay are set to zero, and those for 
OvliB and Qvh00 are smaller, with yp factors of order 
^ RIO for all nuclei. 


and compare them with the unquenched values with gA = 

1.269 also given in Table IXVTl Quenching of gA appears 
to have a major effect on the calculated half-lives, by 
multiplying them by a factor of 6-50.The question of what 
value of gA one needs to use is thus a major concern 

which needs to be addressed. This concern has also been ACKNOWLEDGMENTS 

discussed recently in Refs. (iM 1^. 


VIII. CONCLUSIONS 

In this paper, we have introduced a method for restora¬ 
tion of the isospin properties of the Fermi transition op- 
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of Energy (Grant No. DE-FG-02-9IER-40608), Chilean 
Ministry of Education Fondecyt (Grant No. 1120462), 
and Academy of Einland (Project 266437). 
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Mass number 


FIG. 8. (Color online) Expected half-lives for {m^,) = 1 eV, 
qa = 1.269 and IBM-2 isospin restored nuclear matrix ele¬ 
ments. The figure is in semilogarithmic scale. 


TABLE XVI. Predicted half-lives in OuPP decay with unquenched and maximally quenched gA, 2 and obtained 

from 2vPj3 decay and {m^} = leV. 


Decay 



~0i/ 

U/2 

(lO'V) 




IBM-2 



ISM 


unquenched 

maximally quenched 

unquenched 

maximally quenched 


1.33 


21.5 

13.9 


89.2 

76Ge^76ge 

1.95 


44.0 

8.65 


69.1 

®2Se-)>®2Kr 

0.71 


16.9 

2.22 


18.5 

96Zr^96Mo 

0.61 


16.3 





0.36 


9.8 




iiopd_^iiocd 

1.27 


37.6 





0.63 


19.2 




124Sj,^124pg 

1.09 


35.2 

2.73 


27.6 

i28Te->i28xe 

10.19 


335 

33.5 


344.4 

i30pg^i30xe 

0.52 


17.2 

1.70 


17.6 


10.23 


348 





0.74 


25.5 

2.39 


25.3 

148Nd^l48Sni 

1.87 


68.2 




isoNd^isogm 

0.22 


8.3 




154Sm^454Qd 

4.19 


158 




i60Gd^4®0Dy 

0.63 


24.4 




198pt^l98jjg 

2.77 


125 




232Th_ 

0.44 


22.4 




238u^238p^ 

0.13 


6.7 
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